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The Evolution of type 1 AGN in the IR (15/im). 
The view from ELAIS-S1 * 
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ABSTRACT 

We present the luminosity function of AGN1 (QSO + Seyfert 1). Our sample 

of 21 high-redshift sources is selected from the Preliminary Analysis catalogue in the 
SI field of the European Large Area ISO Survey (ELAIS). To study the cosmological 
evolution of the AGN1 luminosity function, our sample has been combined with a 
local sample of 41 sources observed by IRAS. We find that the luminosity function 
of AGN1 at 15/^m is fairly well represented by a double-power-law-function. There 
is evidence for significant cosmological evolution consistent with a pure luminosity 
evolution (PLE) model i(z)oc(l + z) kL , with kr, = 3.0 — 3.3. The value of kr, depends 
on the existence or not of an evolution cut-off at redshift ~2, and on the adopted 
cosmology. From the luminosity function and its evolution we estimate a contribution 
of AGN1 to the Cosmic Infrared Background (CIRB) of vl v ~ 6 x Kr n W m~ 2 sr" 1 
at 15/im. This corresponds to ~ 2 — 3% of the total observed CIRB at this wavelength. 
Under the usual assumptions of unified models for AGN, the expected contribution of 
the whole AGN population to the CIRB at 15/xm is 10 - 15%. 
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1 INTRODUCTION 

The measurement of the luminosity function (LF) of Active 
Galactic Nuclei (AGN) and its evolution provides funda- 
mental information on the accretion history of the Universe, 
the physics related to it, how structures have been formed 
and on the contribution of this population to the Cosmic 
Background. 

Historically, AGN have been classified using their op- 
tical characteristics and divided into two categories: type 1 
(AGN1) and 2 (AGN2), according to the presence or absence 



* Based on observations with the Infrared Space Observatory 
(ISO). ISO is an ESA project with instruments funded by ESA 
Member States (especially the PI countries: France, Germany, 
The Netherlands and the United Kingdom) and with the partic- 
ipation of ISAS and NASA. 



of broad emission lines in their optical spectra (we will use 
this definition throughout this paper). 

Unified models assume a dusty, IR emitting, torus 
around the central engine. These models distinguish AGN1 
and AGN2 depending on the torus orientation; AGN1 are 
galaxies with low absorption from the mid-IR to the soft 
X-ray while AGN2 are objects with the largest column den- 
sities (e.g. Antonucci 1993; Urry & Padovani 1995). 

The strength of the AGN1 LF evolution has been de- 
termined in the optical, soft and hard X-rays (e.g. Boyle et 
al. 2000; Miyaji et al. 2000; La Franca et al. 2002), while in 
the mid infrared (5-30 /jm: mid-IR) it is still unknown. The 
only measure available at these wavelengths, for AGN1, is 
the local LF at 12 /jm from IRAS data (Rush, Malkan & 
Spinoglio, 1993; RMS hereafter). 

The situation for AGN2 is quite different because they 
are strongly absorbed and thus difficult to select and/or 
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spectroscopically identify. Statistically significant samples of 
AGN2 need to be built up in the far infrared or the hard X- 
rays where absorption is less significant. However, in these 
cases the optical spectroscopic identification is also a diffi- 
cult task. 

For these reasons the contribution of AGN to the Cos- 
mic Infrared Background (CIRB) is still unknown. The only 
estimate is based on AGN1 evolution measured at other 
wavelengths, under assumptions on the AGN spectral en- 
ergy distribution (SED), and the AGN1/AGN2 ratio (e.g. 
Granato, Danese & Franceschini 1997). 

We have spectroscopically identified 25 AGN1 at 15pm 
from the Preliminary Analysis (PA) catalogue of the Euro- 
pean Large Area ISO Survey (ELAIS). Our spectroscopic 
sample is the only one able to probe the high redshift popu- 
lation of AGN1 in this waveband. Other published ISO sam- 
ples (i.e. CFRS, Flores et al. 1999; ISO-HDF north, Aussel 
et al. 1999; ISO-HDF south, Oliver et al. 2002) do not have 
enough sources and/or spectroscopic identifications. 

Here we present our measurement of the evolution of LF 
of AGN1 at 15pm using our data combined with the 12pm 
local sample from RMS. The data are presented in section 
2. In section 3 we describe the method used to derive the 
LF, while the results are presented and discussed in section 



n =0.3 n. =0.7 



2 THE SAMPLE 

ELAIS is the largest single open time project conducted by 
ISO (Oliver et al. 2000), mapping an area of 12 deg 2 at 
15pm with ISOCAM and at 90pm with ISOPHOT. Four 
main fields were chosen in order to reduce the effects of 
cosmic variance (Nl, N2, N3 in the north hemisphere and 
SI in the south) at high Ecliptic latitudes (|/3| > 40°) on 
selected areas of low extinction. 

A Preliminary Analysis (PA) catalogue was produced 
for the 4 main fields (Serjeant et al. 2000). The Final Anal- 
ysis (FA) uses the Lari technique and the catalogue of the 
SI field (J2000,a : 00 h 34"*44.4 s , 5 : -43°28'12", covering 
an area of 3.96 deg 2 ) has been recently released by Lari et 
al. (2001). 

Optical identifications were possible due to an exten- 
sive R-band CCD survey, performed with the ESO/Danish 
1.5m and the ESO/MPE 2.2m telescopes. The spectro- 
scopic follow-up program was carried out with the AAT 
at the Anglo- Australian Observatory (AAO), the 3.6m and 
NTT at ESO/La Silla. The selection of targets for spec- 
troscopic identification was carried out randomly in order 
to cover uniformly ~70% of the sources in the magnitude 
range 17<R<20. The dominant class is star-forming galax- 
ies (about 50%), however AGN (type 1 + 2) constitute a 
significant fraction of the identifications in this magnitude 
range (~30%). The ELAIS-S1 field has also been completely 
covered in the radio at 1.4 GHz down to 0.3 mjy (Gruppi- 
oni et al. 1999), and in 50% of its area in the X-rays with 
BeppoSAX (Alexander et al. 2001). 

In order to assess our sample of optically identified PA 
sources, we took advantage of the already completed FA 
catalogue in the ELAIS-S1 region (Lari et al. 2001). The FA 
catalogue contains 462 sources at 5a confidence level, while 
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Figure 1. Distribution on the Luminosity-redshift plane of the 
type 1 AGN sources used in computing the luminosity function. 



the PA catalogue we used contains 762 sources, as it includes 
also less significant sources. 

All the PA sources with spectroscopic identifications 
were analyzed, and the flux recomputed with the Lari Tech- 
nique at > 4cr confidence level: 23 out of 25 PA AGN1 were 
confirmed and used for our statistical analysis. 

The fluxes of the PA catalogue were recomputed accord- 
ing to the FA calibration by using all the objects in common, 
plus the PA objects with existing spectroscopy and detected 
also with the Lari Technique at > 4<r confidence level. The 
effective a function of flux, of our PA selected sample 

has been statistically estimated by comparing the PA "cor- 
rected" counts (i.e. including only those detected at > 4<j 
confidence level with the Lari Technique) with the counts 
derived from the FA catalogue by Gruppioni et al. (2002). 

Our final sample consists of 21 AGN1 (0.3<z<3.2, see 
Figure 1) with /i5 Mm >lmJy and 17<R<20 (two objects 
fainter than R=-20 fall outside these limits). 

To represent the local mid-infrared population of type 1 
AGN we have combined our data with a sample of 41 AGN1 
belonging to the RMS sample selected from the IRAS Faint 
Source Catalog, Version 2 (Moshir et al. 1991). This sample 
is complete down to 300 mjy at 12pm. 



3 THE LUMINOSITY FUNCTION 

In this work we define L15 as vL v at 15pm. Detailed k- 
corrections (over the 15pm LW3 filter from ISO and the 
12pm filter from IRAS) are necessary in order to derive 
the 15pm luminosities for our sources, as well as to convert 
the 12pm luminosity from the RMS sample into a 15pm 
luminosity. We adopted the mean spectral energy distribu- 
tion (SED) of radio-quiet QSOs from Elvis et al. (1994) as 
representative of our AGN1 population. The standard k- 
correction was computed following Lang (1980). Figure 1 
illustrates the luminosity-redshift space distribution of all 
AGN1 used in this study. 

A parametric, unbinned maximum likelihood method 
was used to fit simultaneously the evolution and luminosity 
function parameters (Marshall et al. 1983) of the combined 
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Figure 2. The apparent R-magnitudes as a function of 15/im 
flux, according to the SED from Elvis et al. (1994) and R-band 
fc-corrections from Natali et al. (1998). Continuous lines represent 
the central values at redshifts 2=0, 1 and 2. Dashed lines repre- 
sent the maximum spread at lrr confidence level. Vertical and 
horizontal lines are the IR and optical limits of the ELAIS-S1 
survey. Plotted squares correspond to our AGN1. 



sample. Since ELAIS identifications were carried out in the 
R-band magnitude interval 17.0<R<20.0, a factor &(z,L), 
which takes into account the optical limits, was introduced 
into the likelihood function. ®(z,L) represents the proba- 
bility that a source with a given luminosity at 15^im (Lis) 
and redshift z, had a R-magnitude within the limits of the 
sample, and was derived taking into account the la internal 
spread in the assumed SED (Elvis et al. 1994). The R-band 
fc-corrections were taken from Natali et al. (1998). In or- 
der to show the dependence of our sample completeness on 
the assumed SED we plot the expected relation (and its la 
spread) between the IR and optical fluxes (Figure 2). From 
this we estimate that our sample of AGN1 is about 76% 
complete. The missing fraction of objects is expected to be 
found in the interval 20.0<R<22.0. The function to be min- 
imized can therefore be written as S = —2 In j£, being C the 
likelihood function, 



S = 



N 



+ 2 



$(2,L)fi(z,L)0(2,L)^cfedL 
dz 



where N is the total number of sources in the two samples, 
0(2, L) is the fractional area of the sky over which a source 
with luminosity L 15 and redshift z could have been observed, 
(dV/dz) is the differential comoving volume and &(z,L) the 
space and luminosity density of the sources being considered. 

Following studies of AGN1 in the optical and in the X- 
ray, we adopted a smooth double power-law for the space 
density distribution in the local Universe, 



d$(L 15 ,z = 0) 
dlogLi 5 



(Lis/Lis)* + (L 15 /L* 15 y 



and tested a Pure Luminosity Evolution (PLE) model of the 
form: 



L 15 (z) = L 15 (0)(1 - 
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Figure 3. Best PLE luminosity function fit (see also table 1). 
Two redshift intervals are shown. The points correspond to the 
space densities of the observed sources, corrected for evolution 
within the redshifts intervals. Sources with 2=[0.0,0.2] are repre- 
sented by diamonds and correspond to the RMS sources. Sources 
with z=[0.2,2.2] are represented by squares, and represent the 
ELAIS-S1 population. Lines plotted are the mid-LF at the cen- 
tral redshift of the interval considered. The error bars are based 
on Poisson statistics at the 68% confidence level. 



Confidence regions for each parameter have been ob- 
tained by computing AS(= A^ 2 ) at a number of values 
around the best-fitting (S m in) parameter, while allowing 
the other parameters to float (see Lampton, Margon & 
Bowyer 1976). The 68% confidence levels quoted correspond 
to AS = 1.0. The normalization factor $* was determined 
by requiring agreement with the observed total number of 
sources (ELAIS + RMS). 

Bidimensional Kolmogorov-Smirnov (2D-KS) probabil- 
ities (see Peacock 1983; Fasano & Franceschini 1987) were 
derived for each of the best fitting PLE models. 

We computed the LF both in a (n m ,fi A ) = (l. 0,0.0) 
and (fi m ,S7A)=(0.3,0.7) cosmology, assuming Ho = 75 
kms -1 MpcT 1 . 



4 RESULTS AND DISCUSSION 

Figure 3 and Table 1 show the results of the fit for the 2 cos- 
mologies considered together with the 2D-KS probabilities. 
The points plotted in Figure 3 correspond to the space den- 
sities of the observed sources, corrected for evolution within 
the redshift intervals as explained in La Franca & Cristiani 
(1997). 

Our (f2 m ,f2A) = (LO,0.0) local LF is not fully consistent 
with the estimate of RMS. Our bright slope f3 ~ 2.9 is 
steeper than the estimate (3 ~ 2.1 of RMS. This difference 
can probably be attributed both to a steeper high-z LF and 
to the effects of evolution on the brightest bins used by RMS. 

The Mid-Infrared LF, for either cosmologies, is steeper 
than that observed in the optical and in the soft and hard 
X-ray. The faint 15/xm slope is a ~ 1.1 while the bright slope 
is f3 ~ 2.9. In the optical Boyle et al. (2000) found a ~ 0.5 
and (3 ~ 2.5 (here we quote slope-1 in order to convert to 
d«I>/dlog L units). In the soft and hard X-rays the faint slope 
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Table 1. Parameter values of the fit of the Luminosity Function 
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Errors correspond to 68% confidence level for a single interesting parameter. 

a. L* corresponds to uL v and is given in ergs . 

b. The normalization, log**, is given in Mpc -3 logL - x . 

c. Contribution (ul v ) of type 1 AGN to the 15/xm Mid-IR background in units of 10 _11 Wm — 2 sr — 1 . 



is a ~ 0.7 while the bright slope is (5 ~ 2.0 (Miyaji, Hasinger 
& Schmidt 2000; La Franca et al. 2002). 

According to the blue optical spectra of our identified 
AGN1, we can exclude the possibility that a different pop- 
ulation has been selected. As a consequence this steepening 
at longer wavelength of the LF, taken at face value, would 
imply a dependence of the spectral energy distribution of 
AGN1 on luminosity and/or redshift. For example, a larger 
IR contribution from the SED of the fainter AGN1 (possi- 
bly caused by a more relevant contribution from the host- 
ing galaxy) could be the origin of the observed behaviour. 
A more detailed discussion on this issue will be carried out 
when the results from the identifications of the FA catalogue 
from ELAIS are available. 

Similarly to what has been observed in the optical, the 
evolution at 15/im is fairly well represented by a PLE model 
with k = 3.0 — 3.3. Our data are not statistically significant 
enough to probe the presence of a cut-off in the evolution 
(z C ut) as observed in the optical and X-ray bands. As shown 
in Table 1, our fits are consistent with a typical value z cut = 
2.0, similar to what has been found from optical and X-ray 
studies (e.g. Boyle et al. 2000; Miyaji et al. 2000). All the 
fits show a good 2D-KS probability, with larger values in a 
(n m ,n A ) = (0.3,0.7) Universe. 

Our estimates of the 15/zm LF allow us to compute 
the contribution of AGN1 to the 15/^m CIRB. Direct mea- 
surements of the CIRB in the 5-100/^m range are difficult. 
Hauser et al. 2001, report an upper limit of 470 and 500 
nWm _2 sr~ 1 at 12 and 25/im respectively, based on the at- 
tenuation of the 7-ray photons (see also Hauser et al. 1998) . 
On the other hand, a lower limit of 3.3 ± 1.3 nWm~ 2 sr -1 at 
15^im is found from the integrated light of source counts de- 
rived from observations with ISOCAM (Altieri et al. 1999). 

The intensity of the Cosmic Background at 15/xm is 
given by, 

Iv = In J dLl5 J dz^ f v (L 15 ,z)$(L 15 ,z) 

where /„(Lis,z) represents the observed flux density of a 
source with an intrinsic luminosity L15 at redshift 2. Ac- 
cording to our derived PLE models, we estimate a contri- 
bution of the AGN1 to the CIRB at 15/^m (in units of vl v ) 
in the range 5.2 — 6.7 x 10 -11 Wm _2 sr~ 1 , depending on 
the cosmology adopted and on the existence or not of a 



redshift cut-off (z cut , see Table 1). The integral has been 
computed for \ogLis(z = 0) > 42 up to z — 5.0. These val- 
ues correspond to about 2 — 3% of the integrated light mea- 
sured by Altieri et al. (1999). If we assume that the ratio 
of AGN2 to AGN1, which is about 4 locally (e.g. Lawrence 
1991; Maiolino & Rieke, 1995), does not change substan- 
tially with redshift and that the shape and evolution of the 
luminosity function of AGN2 is similar to those of AGN1, 
the total contribution for AGN to the background measured 
by Altieri et al. (1999) would be ~ 10 - 15%. This value is 
in rough agreement with what is predicted by the models of 
dusty tori in AGN from Granato et al. (1997). The predic- 
tions from Xu et al. (2001) on the contribution of AGN to 
the CIRB at 15/^m is about a factor of 2 lower than our esti- 
mate. This difference is probably due to two reasons. First, 
the slope of their local LF of AGN1+AGN2 is flatter than 
our estimates at faint luminosities. Second, their assumed 
evolution (a PLE with /cl=3.5 up to z=1.5 and then drop- 
ping as kL~-3) has a lower z cut and a faster decay at high 
redshift. 
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